Initial fluctuation effects on harmonic flows in high-energy heavy-ion collisions 
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Within the framework of a multi-phase transport model, harmonic flows v n (n = 2, 3 and 4) are 
investigated for Au + Au collisions at ^/snn = 200 GeV and Pb + Pb collisions at ^/sjvjv = 2.76 TeV. 
The event-by-event geometry fluctuations significantly contribute to harmonic flows. Triangular 
flow (v's) originates from initial triangularity (£3) and is developed by partonic interactions. The 
conversion efficiency (v n /e„) decreases with harmonic order and increases with partonic interaction 
cross section. A mass ordering in the low pr region and number of constitute quark scaling in the 
middle pr region seem to work roughly for n-th harmonic flows at both energies. All features of 
harmonic flows show similar qualitative behaviors at RHIC and LHC energies, which implies that 
the formed partonic matters are similar at the two energies. 

PACS numbers: 12.38.Mh, ll.10.Wx, 25.75.Dw 
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I. INTRODUCTION 

Results from the Brookhavcn Rclativistic Heavy-Ion 
Collider (RHIC) indicate that a strongly-interacting par- 
tonic matter has been created in rclativistic nucleus- 
nucleus collisions [l|. A powerful probe exposing the 
characteristics of new matter, elliptic flow, has been mea- 
sured via the second Fourier coefficient (i^) in the az- 
imuthal distribution of final particles. It is translated 
from an early stage coordinate space asymmetry, which 
can reflect how the hot matter evolves hydrodynami- 
cally [Ml. 

The V2 data show remarkable hydrodynam- 
ical behaviors, which implies the formed matter is ther- 
malized in a very short time and expands collectively as 
a perfect-like liquid with a very small shear viscosity over 
entropy density ratio (rj/s) [j] Q. Elliptic flow (?j 2 ) has 
been studied widely as functions of centrality, transverse 
momentum (pt) and pscudorapidity (77) etc. A mass- 
ordering at low pt and a Number of Constituent Quark 
(NCQ) scaling at intermediate px for V2 have been ob- 
served, which suggests that a thermalized partonic mat- 
ter is formed and a collective motion is developed prior 
to hadronization @-[T2|- On the other hand, a geometry 
(participant eccentricity) scaling was observed for V2 fluc- 
tuations, which implies not only participant eccentricity 
is responsible for elliptic how, but also the event-by-event 
initial state geometry fluctuations contribute to harmonic 

flow [IM1- 

It has been recently found that the triangular flow (113) 
is not zero in the azimuthal distribution of final parti- 
cles. In fact, because of the non-smooth profile, coming 
from the event-by-event fluctuations of participant nu- 
clcons, it shows a triangular initial geometry shape can 
be transferred into momentum space by hydrodynamical 
evolution. In recent studies, it has been demonstrated 
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that triangular flow significantly contributes on the near- 
side ridge and away-side double bumps in two-particle 
azimuthal correlations jl6l E3. As a new probe, trian- 
gular flow is believed to provide more information about 
the formed hot and dense matter. It has been studied as 
functions of centrality, transverse momentum, pseudora- 
pidity (77) , as well as the relations with the initial triangu- 
larity (£3) and shear viscosity over entropy density ratio 
[Ra - |2T)I ]. However, the dependence of triangular flow on 
the elastic two-body partonic scattering cross section is 
absent. In addition, a possible NCQ-scaling, which has 
been found held by the elliptic flow [2l[, have not been 
studied in details for other v n (n=3,4...) when the initial 
fluctuations are taken into account . 

This work presents the initial deformation scaling of 
elliptic (1*2), triangular (7J3) and quadrangular flows (1*4) 
for different cross sections within the framework of the 
AMPT model [H, H3|. The mass ordering at low p T 
and constituent quark number scaling at higher px for 
the v n are investigated after considering the event-by- 
event initial state geometry fluctuations at RHIC and 
LHC energies. Meanwhile, a special care is discussed for 
s-quark and <f> meson for 7j„-scaling. 

The paper is organized in the following way. A brief 
description of the AMPT model is introduced in Sec. II. 
The results and discussions are presented in Sec. III. 
Finally, a summary is given in Sec. IV. 



II. BRIEF DESCRIPTION OF AMPT MODEL 

A multi-phase transport (AMPT) model consists of 
four main components: the initial condition, partonic 
interactions, conversion from partonic to hadronic mat- 
ter, and hadronic interactions. The initial condition, 
which includes the spatial and momentum distributions 
of minijet partons and soft string excitations, is ob- 
tained from the Heavy Ion Jet Interaction Generator (HI- 
JING) model. Scatterings among partons are modeled by 
Zhang's Parton Cascade (ZPC) model, which includes 
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only two-body scatterings with cross sections obtained 
from the pQCD calculations with screening mass. In the 
default version of AMPT model, partons only include 
minijet partons, and recombine with their parent strings 
when they stop interactions, then the resulting strings are 
converted to hadrons by using the Lund string fragmen- 
tation mechanism. While in the version with the string 
melting mechanism, partons include minijet partons and 
partons from melted strings. And a quark coalescence 
model is used to combine partons into hadrons. The 
dynamics of the subsequent hadronic matter is then de- 
scribed by a relativistic transport (ART) model. Details 
of the AMPT model can be found in a review [23| . Pre- 
vious AMPT calculations have found that elliptic flow 
can be built by strong parton cascade [HHIfl and jet 
losses energy into partonic medium to excite an away- 
side double-peak structure [TH, H3, [H| . Therefore, par- 
tonic effect can not be neglected and the string melting 
AMPT version is much more appropriate than the default 
version when the energy density is much higher than the 
predicted critical density . In this work, we use the ver- 
sion of AMPT model with the string melting mechanism 
to simulate Au+Au collisions at t^/snn = 200 GcV as 
well as Pb + Pb collisions at ^s^n = 2.76 TeV. Since 
collective flow has been built up after the expansion of 
partonic stage, we neglect the final hadronic rescattering 
effects on harmonic flows in this work. 



III. RESULTS AND DISCUSSIONS 

A. Brief definition of v n with initial fluctuations 

We know harmonic flows are defined as the n-th 
Fourier coefficient v n of the particle distribution with re- 
spect to the reaction plane. However, after considering 
event-b y-ey ent fluctuations in the initial density distri- 
bution [16|, the particle distribution should be written 
as 



dN 

— oc 1 + 2 v n cos[n(cf> -ip n )], 



(1) 



where (p is the momentum azimuthal angle of each 
hadron. ip n is the n-th event plane which varies due to 
event-by-event fluctuations and can be calculated by 



arctan 



(r 2 sin(n<£>)) 
(r 2 cos(nip)} 



(2) 



where r and tp are the coordinate position and azimuthal 
angle of each parton and the average (■ ■ ■ ) is density 
weighted in the initial state, and and the superscript r 
denotes initial coordinate space. The n-th order eccen- 
tricity e n for initial geometric distribution is defined as 



yj (r 2 cos(n(/?)) 2 + (r 2 sui(rup)Y 

V) 



(3) 



There is some arbitrariness in the definition of ipn{r) 
and e„ [1^], because one could, for instance, replace r 2 
with r n in Eq. (2) and (3) [H. With this replacement, 
however, 113 and V4 only change little in our calculations 
(less than 3% for w 3 and 12% for V4, respectively, for 0- 
80% centrality). In the following calculations, we will use 
Eq.([2]) and Eq. (J3J) to decide ip n and e n . 

After tp n is determined, the n-th harmonic flow v n can 
be obtained by 



H4>-Vn)]) 



(4) 



In alternative way, tp„ and v„ can be also calculated in 
momentum space as, 



n 



arctan 



(p T sin(n^)) 



(j>t cos(n<p)) 



and 



< i = (c OS [n(0-^)]) 



(5) 



(6) 



where pt and (p are the transverse momentum and az- 
imuthal angle of each hadron, respectively, which is se- 
lected from pseudorapidity |n| > 1 in the final state to 
avoid autocorrelation, and the superscript p denotes final 
momentum space. 

For determined by the final momentum phase space, 
we can even-by- event correct v n into by 



'«S' (e) -*; ,(e) «in[n(Vfi ,(e) -< ,(e) )]' 



tin 



(7) 



where the superscript (e) denotes "event -wise , Rn ^ — 
cos[n(il)n^ — ipn^)] is event-wise event plane resolu- 
tion, and = sin[n{<f> — ipn )] is event-wise sin- 
term harmonic coefficient. We found that the contri- 
bution from sin-term is only approximately 10% , there- 
fore we neglect the sin-term and correct into v n by 

v' n = (vn ^ /Rn ^ ), which is more operable experimen- 
tally. 

It is essential to check if the v n (pr) calculated by differ- 
ent ip n defined in coordinate space and momentum space 
is similar or not, because the determination of ip n in co- 
ordinate space by the Eq.([5|) is not accessible in experi- 
ment. The pr dependences of V2 and V3 with respect to 
ip n determined by initial coordinate and final momentum 
spaces are shown in Figure [1] together with the PHENIX 
V2 data [3l|. We observed that v-i and V3 determined 
by the final momentum space is very close to the ones 
with respect to initial coordinate space. (Note: we check 
that the differences are due to sin— term contributions in 
event-by-event resolution corrections.) Also, the values 
can basically fit the PHENIX data, especially at low pt- 

Based upon the above observations on v-i and V3 with 
different phase space methods, we conclude that they ba- 
sically can present the same results. Therefore in our fol- 
lowing calculations, wc apply the initial coordinate space 
to calculate i\f n and then obtain the corresponding v n - 
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FIG. 1: (Color online) Comparison for V2 (solid symbols) and 
i>3 (open symbols) calculated by tp2,3 denned in coordinate 
space (circles) or tp^ 3 momentum space (squares) in Au + 
Au collisions (0-80% centrality) at ^/snn = 200 GeV from 
the AMPT simulation (3 mb). Note that vi and w$ defined in 
momentum space (squares) have been even-by-event corrected 
by the event plane resolution. The PHENIX data are shown 
by solid stars [3l| . 



B. Initial fluctuations and ratio of v„/e r , 
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FIG. 2: (Color online) e„ and v n in mid-rapidity as functions 
of impact parameter for Au + Au collisions at ^snn = 200 
GeV from the AMPT simulation with partonic interaction 
cross section of 3 mb. 



The ratio of elliptic flow to eccentricity {1)2/ £2) has 
been found to be sensitive to the freeze-out dynamics, 
the equation of state (EOS) and viscosity, however, 
can give more information [1(1 l32T - [36l | . Figure [D shows 
the initial n-th order eccentricity e n and final harmonic 
flow v n (n = 2, 3 and 4) in mid-rapidity as a function 
of impact parameter for Au + Au collisions at y / sjv"/v = 
200 GeV from the AMPT model simulations. The elastic 
two-body scattering cross section in the parton cascade 
process is set to be 3 mb. 

As presented in Figure [2j the n-th order eccentricity 
e n (n = 2, 3 and 4) increases with impact parameter. £2 
is larger than £3 and £4, except in very central collisions 
where e n looks similar to each other. It is consistent 
with the trend given by Lacey et al. who applied MC- 
Glauber model, but gave a little smaller magnitude for 
peripheral collisions [37] . On the other hand, the coeffi- 
cients of anisotropic flow v n (n — 2, 3 and 4) show rising 
and falling with impact parameter. Also, v n has a larger 
magnitude for lower harmonic than higher harmonic. 

Similarly, Figure [3] shows the initial geometry defor- 
mation £„ (n = 2, 3 and 4) for Pb + Pb collisions at 
y/sNN = 2.76 TeV, which demonstrates very similar be- 
haviors as the RHIC energy. 

Once we have v n and £„, we can discuss the ratio v n /e n 
(n = 2, 3 and 4). Figure 2] shows impact parameter and 
partonic cross section dependences of v n /e n for Au + Au 
collisions at y/sjyw = 200 GeV. The value of ratios de- 
creases with impact parameter, which implies that the 
conversion from the initial geometry asymmetry to fi- 



nal momentum anisotropy is less efficient for peripheral 
collisions than for central collisions. And for higher har- 
monics, there is also less conversion efficiency. The trend 
for v n /e n as a function of impact parameter looks sim- 
ilar for the different partonic interaction cross sections 
of 3, 6 and 10 mb. However, the magnitude of v n /£ n 
decreases with the cross section, which reveals that the 
conversion from the initial geometry asymmetry to the fi- 
nal momentum anisotropy becomes weaker for a smaller 
cross section. This indicates that frequent parton-parton 
collisions help the system to develop the harmonic col- 
lectivity. 



From Figure [H we also saw that the v n /e n becomes 
smaller for higher harmonic order, this may reflect the 
viscous damping. Recently, it was claimed that the rela- 
tive magnitude of the higher-order harmonics (v n ,n > 3) 
can provide additional constraints on both the magni- 
tude oi n/s and the determination of initial condition 
[H], H3, [H| . Fig. [5] shows the n-dependencc of v n /e n for 
Au + Au collisions at ^/Inn = 200 GeV and Pb + Pb 
collisions at v / ijv r /v = 2.76 TeV for 0-20 % centrality and 
low p T region (p T < 0.55 GcV/c) (3 mb) with corre- 
sponding exponential fitting functions. Compared with 
PHENIX data , only the trend can be reproduced. 
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FIG. 5: (Color online) Harmonic order dependence of v n /e n 
for Au + Au collisions at ^snn = 200 GeV (circles) and Pb + 
Pb collisions at v /s^Fiv = 2.76 TeV (squares) from the AMPT 
simulation (3 mb). PHENIX data are shown by triangles. 
The centrality is 0-20% and p T < 0.55 GeV/c. The curves 
are exponential fitting functions. 



FIG. 3: (Color online) Same as Fig. but for Pb + Pb colli- 
sions at ^/sjvjv =2.76 TeV. 
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FIG. 4: (Color online) v n /e n (circles for n=2, squares for n~3, 
and triangles for n=4) as a function of impact parameter for 
Au + Au collisions at -y/ijviv = 200 GeV in mid-rapidity for 
different parton interaction cross sections (10 mb, 6 mb and 
3 mb). 



C. pr dependence of v n with different partonic 
cross sections and comparisons with the data 

Figure |5] presents our simulations of V2, v$ and v± as 
a function of px with different parton interaction cross 
sections together with the PHENIX data [|||. For tri- 
angular flow, it totally arises from the event-by-event 
fluctuations of the initial collision geometry, because it 
persists zero if without considering the fluctuations. v n 



(n = 2, 3 and 4) decreases when parton-parton cross sec- 
tion decreases. Experimental data of i>2 can be described 
by the large cross sections (from 3 mb to 10 mb), after 
one considers of the initial fluctuations. However, the 
AMPT model underestimates the data if without tak- 
ing the initial fluctuations into account. Recently, Xu 
and Ko adjusted more parameters in the AMPT model, 
which include not only parton interaction cross section 
but also the paramctrization of the Lund string fragmen- 
tation, and found that a smaller cross section of 1.5 mb 
is good to describe both the charged particle multiplicity 
and elliptic flow |40j . In our work, we will not focus on 
how to further improve parameters, but we do find that 
initial geometry fluctuations significantly affect harmonic 
flows and should not be ignored. 

The transverse momentum dependences of i>2 and V3 
with different cross sections in four different centrality 
bins are shown in Figure [7J and M The PHENIX data 
is also accompanied [4l|. For each centrality bin, v-2 
and increase with the cross section. For elliptic flow 
(Fig. [7]), data seem to prefer a bigger cross section in 
higher transverse momentum range. In the case of tri- 
angular flow similar trend is present in Fig. [8l though 
V3 shows a less centrality dependence than V2, which is 
consistent with the trends shown in Figure 

The transverse momentum dependences of V2 and V3 
are also calculated for four different centrality bins in Pb 
+ Pb collisions at ^s^n = 2.76 TeV for LHC energy, 
which are shown in Figure [9] and [10] together with the 
ATLAS data [42j . Similarly, V2 and i>3 increase with the 
cross section from 3 mb to 10 mb, which can basically 
describe the ATLAS data, 

Even though a general behavior of p^-dcpcndcnt v n 
can be nicely demonstrated by the comparison of our cal- 
culations with the data, we found that the AMPT Simula- 
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FIG. 6: (Color online) (a)-(c): Transverse momentum depen- 
dences of V2, V3, and V4 for Au + Au collisions (0-80% central- 
ity) at y/SNN = 200 GeV in mid-rapidity for different cross 
sections [10 mb (circle), 6 mb (square), and 3 mb(triangle)], 
where solid symbols are v n with considering initial fluctua- 
tions and open symbols are those without considering initial 
fluctuations. The PHENIX data are shown by solid stars [3||. 



tions can only describe the trend of the data qualitatively. 
Actually, AMPT can not describe the V2 and 113 data si- 
multaneously with a same cross section. For example, 
the 3mb results describe the V2 and V3 for 0-10% central- 
ity but undcrprcdict other centralities at RHIC energy. 
The lOmb data describe the V2 for 10-20% centrality but 
are unable to reproduce w 3 for the same centrality. Also, 
the lOmb results describe the V2 and W3 for 30-40%, but 
underpredict the high pr data for 50-60% centrality. For 
LHC data, no AMPT calculations can describe the t>2 
and V3 for the 0-10% centrality range. 



D. NCQ-scaling of v n at RHIC energy 
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FIG. 7: (Color online) V2 as a function of px in Au + Au 
collisions at v / ijviv" = 200 GeV for different centrality bins 
(0-10%, 10-20%, 30-40% and 50-60%) with different cross sec- 
tions. Circles, squares and triangles represent the calculation 
with 10 mb, 6 mb and 3 mb, respectively. The PHENIX data 
shown by solid triangle [4lj [. 
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FIG. 8: (Color online) Same as Fig. [7] but for V3. 



(v2/n q vs KEt/Uq) reveals a universal scaling of V2 for 
all identified particles over the full transverse kinetic en- 
ergy (KEt) range, which is more pronounced rather than 
Pt @, Iii-ll2l E3| ■ Such scaling indicates that the collec- 
tive elliptic flow has been developed during the partonic 
stage and the effective constituent quark degree of free- 
dom plays an important role in hadronization process. 

For higher even-order harmonics, V4 and vq etc, appear 

1/2 



For elliptic flow, a mass ordering (the heavier the 
hadron mass, the smaller the V2) and a NCQ-scaling 
(baryon versus meson) have been observed at low and 
intermediate px, respectively, in Au + Au collisions at 

y/sNN — 200 GeV [i^. The observed NCQ-scaling however, seems to be scaled by n} q [2l|. It is interesting 



to be scaled as v n oc v 2 [43 and their NCQ-scaling has 
also been suggested in Ref. [461 ] . Even in very low energy 
heavy ion collisions, the ^-scaling and the tti/i>f -scaling 
have been suggested for light nuclear clusters in nucle- 
onic level interaction (4?J . Instead scaling by the number 
of constituent quarks (n„) for v 2 -, the measured data V4, 
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FIG. 9: (Color online) V2 as a function of pt in Pb + Pb 
collisions at -y/5jvjv = 2.76 TeV from the AMPT simulations 
(3 mb and 10 mb) for different centrality bins (0-5%, 10-20%, 
30-40%, and 40-50%) at midrapidity. Circles and triangles 
represent the calculations with 10 mb and 3 mb, respectively. 
The ATLAS data shown by solid triangles [12]. 
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FIG. 10: (Color online) Same as Fig. H but for 



V:i. 



to check if these scaling relations are still valid for the 
v n calculations in which the initial fluctuations are taken 
into account, including the odd harmonics, such as V3. 

Figure [Til presents V2, «3 and of different types of 
hadrons in mid-rapidity for Au + Au collisions (0-80% 
centrality) at ^/snn = 200 GeV. Figure Qj] (a) shows 
that V2 preserves an obvious mass ordering in relatively 
low px region, and hadron type grouping in intermediate 
Pt region, even after considering cvent-by-event fluctua- 
tions. Similarly, and 174, [Figure QT] (b) and (c)] also 
present a mass ordering in the low pt region. The study 
on v n of different hadron species will give more informa- 
tion about the initial geometry and the viscosity of hot 
and dense matter [34[ ■ 

As shown in Figure [T^ a), V2 scaled by the number of 
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FIG. 11: (Color online) (a)-(c): Transverse momentum de- 
pendences of V2, V3, and V4 for different hadron species in Au 
+ Au collisions (0-80% centrality) at y / ijvjv = 200 GeV in 
mid-rapidity from the AMPT simulations (3 mb), with con- 
sidering of initial fluctuations. 



constituent quarks (v2/n q ) as a function of the transverse 
kinetic energy (KEt = \/ p\ + m 2 — to) scaled by the 
number of constituent quarks (K Ex / n q ) shows a univer- 
sal scaling regardless of the initial geometry fluctuations 
are taken into consideration or not. The only difference 
is that the initial fluctuations enhance the value of V2/n q . 
Therefore, the initial fluctuations have little effect on the 
breaking of the NCQ-scaling for elliptic flow. Figure [T^] 
(b) and (c) display v^/n^J 2 and v^/n 2 , for all hadrons 
as a function of KEx/n q , respectively, when the initial 
fluctuations are considered. Form the above results, it 
seems that v n can still be roughly scaled by n q ^ 2 for all 
hadrons as a function of KEr/nq. Of course, the scaling 
behavior is not perfect within the present statistics. For 
example, the amount of spread between different particle 
species is less than 10% for the W2-scaling, it is less than 
20% for the i^-scaling, but it can reach 20-30% for the 
^-scaling. 

In order to understand possible origin of the NCQ- 
scaling of v n for different mesons and baryons, we also 
check v n of u. d and s-quarks as a function of pr or 
KEt- As expected, there exists similar NCQ-scaling of 
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FIG. 12: (Color online) (a)-(c): V2/n q , vz/rh/ 2 and v^/n^ 
as a function of KEr/riq in Au + Au collisions (0-80% cen- 
trality) at ^/F^T = 200 GeV from the AMPT simulations (3 
mb), where the solid symbols are the results for considering 
initial fluctuations (w/), while the open ones are for without 
considering initial fluctuations (w/o). 



v n (n=2-4) for all those constituent quarks. Furthermore, 
we find that the values of v n /n q n ^ 2 of different hadrons 
are similar to the values of v n of u, d, s-quarks, which 
reflects that the NCQ-scaling of v n for different hadrons 
stems from partonic level. 

Furthermore, the ratios of v^/v^ 2 and v^jv\ as func- 
tions of px for three different centrality bins (10-20%, 20- 
30%, and 30-40%) in Au + Au collisions at ^/IWn = 200 
GeV are shown in Figure[T2] It shows that the both ratios 
of v 3 /vl /2 [Figure [H (a)] and w 4 /u| [Figure [U (b)] ex- 
hibit centrality dependences, i.e. more central collisions 
result in more larger ratios. However, it is almost inde- 
pendent of pt for each centrality bin, which is consistent 
with the scaling of v n /riq^ 2 observed in Figure [T2l How- 
ever, it is difficult to obtain the relationship between v n ^ q 
and U2,ij, such as between v 3 and V2, in terms of a sim- 
ple coalescence model in Ref. [48[, because ^3 is purely 
determined by initial geometry fluctuations which is inde- 
pendent of ip2 ■ It is interesting that recently Lacey et al. 
linked such a scaling to the acoustic nature of anisotropic 
flow to constrain initial conditions, rj/s and viscous hori- 



zon [49|. It further gives more insights on the dynamics 
of strongly-interacting partonic matter and constituent 
quark degree of freedom in hadronization process. 



(a) 



• 10-20% 
■ 20-30% 
A 30-40% 



ft|<1.0, 3 mb 



(b) 



CN CM 
> 



0.5 1 1.5 0.5 1 

p T (Gev/c) 



1.5 



FIG. 13: (Color online) The ratios of v 3 /v% /2 (a) and v 4 /d| (b) 
as a function of pr for three different centrality bins (10-20%, 
20-30%, and 30-40%) in Au + Au collisions at ^S^F = 200 
GeV from the AMPT simulations (3 mb). 



E. NCQ-scaling of v n at LHC energy 

At the same time, the mass ordering in the low pt re- 
gion and the NCQ-scaling in intermediate px region of 
v n are also investigated at LHC energy in AMPT sim- 
ulations. Fig. [14] presents the results of V2, V3, and V4 
for different hadron species in Pb + Pb collisions (0-80% 
centrality) at ^/snn = 2.76 TeV in mid-rapidity from 
the AMPT calculations (3 mb), with considering initial 
fluctuations. It displays that the mass ordering is satis- 
fied, i.e. v n decreases from n, k, p, <j> to A in the lower 
Pt region (note that <J) is very close to p in the figure. 
However, the strict mass-ordering needs c/Vs v n is little 
less than p J s v n ). The baryon- meson typing is also ev- 
ident above pr ~ 1.2 GeV. By transformation of pr to 
KEt/Uq as well as v n to v n /n^ 2 , the results of Vi/riq, 
v-ilr?q 2 and v^/n 2 as a function of KEx/nq are shown 
in Fig. [151 Again, the NCQ-scaling of v n is roughly kept 
except for <p meson whose v n is a little larger. Of course, 
the amount of spread between different particle species 
for u„-scaling keeps similar as RHIC energy. Comparing 
with the above v n results at RHIC, LHC results are very 
similar but they reveal larger v n values than RHIC's due 
to stronger partonic interactions at higher energy. But 
in general, the partonic matter formed at LHC energy is 
very similar to that created at RHIC energy. 

As mentioned above, </j-meson shows a little larger 
value of v n /n q n / 2 as compared to those of other hadrons 
in Fig. 1151 Keep in mind that <fr- meson is always a very in- 
teresting hadron in previous studies since its mass is close 
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FIG. 14: (Color online) (a)-(c): The transverse momentum 
dependences of V2, Vz, and V4 for different hadron species in 
Pb + Pb collisions (0-80% centrality) at ^s]7n = 2.76 TeV 
in mid-rapidity from the AMPT simulations (3 mb), with 
considering initial fluctuations. 



to proton but it is a multi-strange meson 0, l50rl54j . This 
could be understood from the parton's v n in the same 
condition: v n of s-quark displays a slight deviation from 
the u (d)-quarks (not shown here). The reason could be 
that the v n of heavier strange quarks has a smaller value 
at low px but a larger value at high px, i.e. the mass 
ordering of partonic flow. However, a larger collective 
radial flow at LHC energy could push heavier s-quark to 
have stronger v n . The effect is of course more distinct at 
LHC energy because of larger initial partonic pressure. 
In contrast, in low energy RHIC run, such as 11.5 GeV/c 
Au + Au collision, s-quark may not reach full thermal- 
ization and therefore result in a less V2 of s-quark as 
compared to u(d)-quarks, which can lead to a smaller V2 
of (f>, i.e. the violation of the ^-scaling for the </>-mcsons 
relative to other hadrons as observed in the STAR data 
[54] as well as in a simulation [ijj ■ Considering that the 
</>-meson is coalesced by ss in the present AMPT model 
calculation, it will certainly induce a larger v n /n q n ^ 2 for 
(j) in comparison with other hadrons, as shown in Fig. 1151 
Unfortunately, the data of eft's v n is not available yet at 
LHC energy, which is worth waiting for checking. 



FIG. 15: (Color online) (a)-(c): V2/n q , vz/rh/ 2 and v^/rig as a 
function of KEr/riq in Pb + Pb collisions (0-80% centrality) 
at ■ s /snn = 2.76 TeV from the AMPT simulations (3 mb). 

Before closing the discussions on the NCQ-scaling of v n 
in this subsection, we remind that the hadronic rescat- 
tcring process is not yet taken into account in our cal- 
culation. Recently, the ALICE data shows that proton's 
V2 and V3 seem to deviate from the NCQ-scaling of v n 
of charged n and K [55[ . The reason could be stronger 
final-state interaction for protons. Detailed model inves- 
tigations are underway. 



IV. SUMMARY 

Within the framework of a multi-phase transport 
model, we investigated the different orders of harmonic 
flows, namely elliptic flow, triangular flow and quadran- 
gular flow for Au + Au collisions at ^Jsnn = 200 GeV 
as well as Pb + Pb collisions at ^/s^vjv = 2.76 TeV when 
the initial geometry fluctuations are taken into account. 
Basically, the harmonic flow is converted from initial ge- 
ometry shape via parton cascade process, and its conver- 
sion efficiency (v n /e n ) decreases with the increasing of 
harmonic order as well as the decreasing of the partonic 
cross section at both RHIC and LHC energies. Depen- 
dences of transverse momentum, centrality and partonic 
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cross section of the v n (n=2, 3 and 4) have been stud- 
ied and compared with data. For each centrality bin, v 2 
and V3 increases with cross section, especially at higher 
transverse momentum. 

Triangular and quadrangular flows also roughly 
present a mass ordering in low px region and the num- 
ber of constitute quark scaling in intermediate pr region, 
similar to the behaviors of elliptic flow. Form our re- 
sults, a NCQ-scaling of v n /n^ 2 versus KE T /n q for dif- 
ferent hadrons holds for harmonic flow (v n , n = 2, 3 and 
4), which can be related to w n -scaling in partonic level. 
From all above results, it implies that the formed par- 
tonic matter should be very similar for RHIC and LHC 
energies. 
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